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^ ■ area except for an extremal limit when this property is preserved 

only for orthogonal charges. Extremal dyon configurations with 
orthogonal charges have finite horizon radii continuously varying 
from zero to the ADM mass. New general solution is endowed 
with a NUT parameter, asymptotic values of dilaton and axion, 
and a gauge parameter which can be used to ascribe any given 
value to one of scalar charges. 
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1 Introduction 



Recently there has been a substantial progress in understanding the statis- 
tical origin of the black hole entropy in the string theory (for a review and 
references see f^j). In achieving this goal it was important to realize that 
black holes in theories with dilaton and multiple vector fields may have BPS 
saturated states with the finite area of the horizon. For toroidally compact- 
ified heterotic string the effective four-dimensional theory is iV = 4 super- 
gravity coupled to Abelian vector multiplets. The above feature is manifest 
already at the level of the pure N = 4 supergravity. Bosonic sector of this 
theory consists of the gravity coupled system of six £7(1) vector fields inter- 
acting with dilaton and axion. It is well known that, with only one U(l) 
component excited, the dilaton causes the horizon of the extremal black hole 
to shrink to the singularity 0, |3], |J. With more than one U(l) charges 
present, some BPS configurations may possess the finite horizon area like 
the Reissner-Nordstrom solution. This was demonstrated by Kallosh et al. 
U within the two-vectors model without axion. Extremal black holes with 
finite horizon area have N — 1 residual supersymmetry H, H] contrary to an 
extremal "dilatonic" black hole with vanishing horizon area, which exhibits 
two supersymmetries unbroken @. Similar conclusions were made in the 
heterotic string effective theory with additional vector multiplets present || . 

In view of the above it is important to investigate the space of classical 
solutions to pure and matter coupled N = 4, D = 4 supergravity in more 
detail. Although many particular black hole solutions were obtained earlier 
using various solution generating techniques, or solving Bogomol'nyi type 
equations in the BPS limit ||, some recent study |J revealed that this the- 
ory still requires further analysis. Of particular importance are stationary 
solutions which can be explored using powerful integration methods devel- 
oped earlier in the vacuum General Relativity and Einstein-Maxwell theory 
|T0|j . Their generalization to dilaton-axion gravity interacting with one vec- 
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tor field was elaborated recently [|11], [12| . The approach consist in derivation 
of the three-dimensional a-model, revealing the isometries of the target man- 
ifold (potential space), and utilizing its Kahler structure to achieve a more 
concise formulation. This procedure has direct similarity with the Ernst ap- 
proach to General Relativity [[R| based on the introduction of complex po- 
tentials parameterizing Kahler target manifolds SL(2,R)/SO(2) (vacuum) 
and SU(2,1)/(SU(2) x U(l)) (electrovacuum) @. It brings considerable 
simplifications into solution generating techniques both at the level of the 
finite-dimensional symmetry group, and the infinite-dimensional Geroch- 
type symmetries if further assumption of the second spacetime symmetry 



(axial) is made [16 



Kahler parameterization of the target manifold of the stationary dilaton- 



axion gravity with multiple vector fields was given in |L7[. The purpose of the 
present paper is to develop the corresponding solution generating technique 
and construct the most general non-rotating black hole solution to iV = 4 
theory in a form manifestly covariant under the three-dimensional ^-duality 
group. We classify isometries of the target manifold in usual terms of Gen- 
eral Relativity (Ehlers-Harrison transformations, electric-magnetic duality, 
gauge, scale transformations and SO(p) rotations), and give explicitly the 
corresponding holomorphic maps (p being the number of U(l) vector fields). 
Choosing particular combinations which preserve asymptotic conditions, we 
apply them to the Schwarzschild solution to covariantize it with respect to 
the [/-duality group SO (2, 2 + p). Thus new solution depending on 2p + 5 
real parameters is generated. 

Our solution includes many particular previously known black hole config- 
urations, and allows one to determine their position in the general parameter 
space. The metric depends on five parameters (including NUT), while the 
material fields are determined by two charge vectors, asymptotic values of 
dilaton and axion, and a gauge parameter which can be used to generate 
solutions with a prescribed value of one of scalar charges. The metric is pre- 
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sented in the curvature coordinates facilitating its physical interpretation and 
analysis of the BPS limit. It is shown that generic non-extremal black hole 
possesses two non-singular horizons and has a modified Reissner-Nordstrom 
form unless charges are parallel. If they are parallel, the internal horizon 
shrinks to the singularity and the metric reduces to that of the "dilaton" 
black hole §]• In the case of orthogonal charges of equal norm the met- 
ric is pure Reissner-Nordstrom (thus we prove that this configuration found 
previously in || is unique). In the general case the metric deviates from 
the usual Reissner-Nordstrom one by a "variable mass" term, asymptoti- 
cally equal to ADM mass and depending on one additional parameter. In 
the BPS limit the (NUT-less) solution with non-orthogonal charges has a 
"dilatonic" form (event horizon shrinking to the singularity). When charges 
are orthogonal but have non-equal norms, the radius of the event horizon in 
curvature coordinates varies between zero value (for the non-dyon case) and 
the ADM mass (for the case of equal norms). 

The plan of the paper is as follows. In the sec. 2 the four-dimensional 
theory is reviewed and the symmetries of the equations of motions are briefly 
discussed. Then we perform 1 + 3 decomposition of the metric and derive a 
cr-model representation in three dimensions (sec. 3). Infinitesimal isometries 
of the target manifold are presented in the sec. 4. Complex parameterization 
of the target manifold is introduced in the sec. 5 where the holomorphic maps 
corresponding to the target space isometries are constructed. The sec. 6 is 
devoted to the covariantization of the Schwarzschild solution with respect to 
[/-duality of the three-dimensional theory. In the sec. 7 we discuss the choice 
of coordinates, analyze various particular black hole configurations and study 
the BPS limit. We conclude with brief remarks clarifying the relationship to 
previously done works. 
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2 Preliminaries 



We consider the gravity coupled system on p Abelian vector fields A a , a = 
1, ...,p, one scalar <fi (dilaton), and one pseudoscalar k (axion) described by 
the action 

S = lL / {~ R + 2d ^ d ^ + - e-^F^F"^ 

-kF^F^} ^9d 4 x, (2.1) 

where F a ^ u = \E^ vXt F% t , F a = dA a , R is the (four-dimensional) scalar 
curvature corresponding to metric g^ v , sum over repeated a is understood. 
For p = 6 this is the bosonic part of the N = 4, D = 4 supergravity. 

Equations for vector fields following from this action have a form of mod- 
ified Maxwell equations 

V v {e~ 24> F a ^ + KF a ^) = 0, (2.2) 

together with the Bianchi identities 

V v F a » v = 0. (2.3) 

The dilaton satisfies the curved space D'Alembert equation with axion 
and vector sources: 

V M V M = - e -^F^F a ^ u + -e 4 *(V/t) 2 , (2.4) 

while the axion is generated by the pseudoscalar invariant of vector fields 

V M (e 4 ^V M /t) = -F^F a " u . (2.5) 

The system is closed by the Einstein equations with vector and scalar sources 

R, v = 2V /1 0V„0 + Vv, K V^ + e-^(2F- x F a \ + ^ T F aXT g, u ), (2.6) 
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where (V/t) 2 = g^ v (y | t /t)(V and is the covariant derivative with 
respect to the 4-dimensional metric g^ v . 

The action (2.1) is invariant under the global SO(p) rotations of vector 
fields, 

ai - n* b A* n a b n b c = 5 a c . (2.7) 



Another symmetry of the equations of motion, known as S-duality [18 
be concisely expressed using the complex axidilaton variable 



can 



z = K + ie 2<t> . (2.8) 

In terms of z the action (2.1) can be rewritten as 

S = -— J {R + 2VzVz(z - z)- 2 - 2Re (izF^J* 1 *) } \^gd 4 x, (2.9) 

where T a = (F a + iF a ) /2, F afJ,v = |e M ^ Ar F^ T (bar denotes complex conjuga- 
tion). Then the (four-dimensional) S— duality transformations read 

z ^az+jl_ {2 + 5) ^ t a5 -f3j^l. (2.10) 

72 + * 

They leave invariant the kinetic term of the 2-field in the action (2.9) as well 
as the full set of equations of motion (but not the action). 

It is worth comparing the present theory with the Einstein-Maxwell one 
(bosonic sector of N = 2 supergravity) : 

s = ~J {R + 2 (^V^)} V^gd'x. (2.11) 

In this case both the Maxwell equations and Bianchi identities can be com- 
bined into one complex equation 

V^ a>MU = 0, (2.12) 

while the Einstein equations read 

R, v = tf%J$. (2.13) 
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Electric-magnetic duality now is an Abelian symmetry of the full action: 

^x - e*\^, (2.14) 

where a is a constant parameter. Notice that the action of the dilaton-axion 
theory (2.9), which has two additional field variables with respect to (2.11), 
possesses at the same time larger symmetries (three-parametric SL(2, R) in- 
stead of one-parametric U(l)). Due to this enhancement of four-dimensional 
symmetries, further three-dimensional reduction of the dilaton-axion gravity 
is describable as a non-linear a-model on a symmetric space similarly to the 
Einstein-Maxwell theory. 

3 3+1 decomposition 

Let us consider the system of equations (2.2-6) on a class of metrics depending 
effectively on three coordinates. In coordinate-independent way this can be 
expressed as a restriction on spacetime to admit a Killing vector field. It 
turns out that, if the Killing field is non-null everywhere, the full system of 
equations can be presented in the form of a gravity coupled non-linear sigma- 
model on the symmetric space 5*0(2,2 +p)/(SO(2) x S0{2 + p)) or one of 
its non-compact form. Aiming to investigate black holes, we assume here 
that the Killing vector is timelike in an essential region of spacetime. Then 
ten components of g^ v can be expressed using the standard Kaluza-Klein 
ansatz through six components of the three-space metric hij, one (co) vector 
u>i, = 1,2,3), and one scalar /, depending on the space coordinates x % : 

ds 2 = g iiV dx ii dx v = f(dt - ujidx 1 ) 2 - jh ij dx i dx j . (3.1) 

To decompose vector fields, one has first to introduce the set of electric 
potentials Aq = v a j\[2 (the coefficient helps to avoid undesired numerical 
factors in the three-dimensional formulation). The spatial parts of four- 
potentials can be traded for another set of scalars, magnetic potentials u a , 
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using the equations (2.2) for v = i. Since spatial part of the tensor at the 
left hand side is divergenceless, one can write 

e -20 F aij + K paij = f e ijk d a (3^) 

V2h 



This relation can be put into an alternative form 

FSl = ^f»l w a k = v(d k u a -Kd k v a ), (3.3) 

and w k can be considered as covariant components of the three-dimensional 
vector w a . 

The remaining Maxwell equations and Bianchi identities now lead to the 
following set of three-dimensional equations for v a , u a : 

v (rv 2 ^ w) + r 2 TVu a - v (rVe 2 ^w a ) = o, (3.4) 

V (/"'e^wj - f~ 2 TVv a = 0, (3.5) 
where three-dimensional vector r dual to two-form duj is introduced 

ijk 

t* = -f-^=d 3 u k . (3.6) 

Here (and in what follows) Latin indices are raised and lowered using the 
three-metric and its inverse h % \ and V denotes a three-space covariant 
derivative. 

To clarify the role of r, the mixed components of the Ricci tensor should 
be invoked |T9| : 

* = (37) 

Upon substitution in the source of the corresponding components of the 
Einstein equations, one finds 

S k d k Tj = 2e ijk d jV a d k u a . (3.8) 
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This equation can be integrated by introducing the NUT (twist) potential x- 

n = d iX + v a d lU a - u a diV a . (3.9) 

Multiplying (3.9) by f~ 2 , and taking into account (3.6), one obtains the 
following second order equation for the NUT-potential: 

f(A X + v a Au a - u a Av a ) = 2(Vx + v a Vu a - u a Vv a )Vf, (3.10) 

where 

A = V 2 = -^diiVhti'dj) (3.11) 

is a three-dimensional Laplacian. From R 00 component of the Einstein equa- 
tions one obtains an equation for the metric function / = g 00 in a similar 
form 

fAf - (V/) 2 + r 2 = / (e 2 V a ) 2 + e- 2<t, (Vv a ) 2 ) . (3.12) 
Finally, the dilaton and axion equations (2.4-5) reduce to 

1 / \ e 4 * 

A0 = — (e" 2 *(W) 2 - e 2 V a ) 2 ) + ^-( Vft 2 > ( 3 - 13 ) 

V(e 40 V/t) = je 2(l) w a Vv a , (3.14) 

Now one can check that this system of equations is derivable from the 
following three-dimensional action 

Sm = l{ w 2 ( (v/)2 + (Vx + yaVua " m£W)2 ) 2 W) 2 + 

+ h^(Vn) 2 - j (e 2 V a ) 2 + e" 2 *(W) 2 ) ^Vhd 3 x, (3.15) 

which describes a non-linear cr-model with 2p + 4 scalar fields ip A = (/, \i 
v a , u a , k ) 0), A = 1, 2p + 4. The remaining (spatial) Einstein equations 
reduce to three-dimensional Einstein equations for the metric hif 

Kij = ^j- 2 (f,f, + nr 3 ) + 20,0, + ^e 4 V^j - j (e^X + **<v>f) , 

(3.16) 
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where IZij is three-dimensional Ricci tensor. Comparing with (3.15), one 
can see that the source term is derivable from the same action (3.15) as the 
energy-momentum tensor. Therefore, the full system of equations follows 
from the three-dimensional gravity coupled a-model action 

S a = J (-K + g AB d iV A d jV B h ij ) Vhd 3 x, (3.17) 

where the target space metric can directly read off from (3.15) as follows 



— [df + (d X + v a du a - u a dv a ) 2 } + 2dcf> 2 



+ -e^dn 2 - - {e 2(t> (du a - ndv a ) 2 + e~ 2(t> (dv a ) 2 ) . (3.18) 

The target manifold is a 2p + 4 dimensional pseudoeuclidean space of 
the signature 2p — 4 which encodes the hidden global symmetries of the sta- 
tionary dilaton-axion gravity with p vector fields (p = 1 theory was studied 
earlier in |l.IJ). 111 absence of scalar fields it reduces (up to generalization to 
arbitrary p) to the potential space of Neugebauer and Kramer for the sta- 
tionary Einstein-Maxwell system [2(J. It is worth noting, however, that the 



target space of the Einstein-Maxwell theory does not constitute a subspace 
of the full target manifold (3.18). Contrary to this, the vacuum Einstein 
target manifold, parameterized by /, x ^ s a subspace. This means that there 
is an intrinsic connection between stationary solutions of the vacuum (but 
not electrovacuum) General Relativity and solutions to the present theory. 



4 U-duality 

Global symmetries of the three-dimensional theory are manifest as isometries 
of the target space. As far as we have an explicit expression for the target 
space metric, we can explore its isometries by solving Killing equations 

K A , B + K B , A = 0, (4.1) 
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where covariant derivatives refer to the metric (3.18). General solution to 



(4.1) was given in J21J for a particular case p — 1. Generalization to arbitrary 
p is rather straightforward in the gauge and .S-duality sectors, but is non- 
trivial in the Ehlers-Harrison sector. Here we give the full set of solutions to 
(4.1) for arbitrary p. 

Gauge transformations 

This sector includes 2p + 1 linear constant shifts of electric and magnetic 
potentials and a NUT potential. No physical quantities are changed. Since 
electric and magnetic variables are mixed with NUT in (3.9), the electromag- 
netic gauge transformations include appropriate variations of x 

v a ^v a + e a , x ^ x - u a e\ (4.2) 

u a -»• u a + m a , x~^X + v a m a , (4.3) 

X^X + 9, (4-4) 

where e a , m a and g are real parameters. The corresponding Killing vectors 
are 

K e a = d v a - u a d x , K™ = d ua + v a d x , K 9 = d x . (4.5) 

Scale and SO(p) rotations 

An invariance of (3.18) under rescaling 

/^e 2 7, x^e 2s X , v a ^e s v a } u a -> e s u a (4.6) 

is obvious, with k and <ft unchanged. This transformation changes the four- 
dimensional metric. The corresponding Killing vector is 

K s = 2fd f + 2 x d x + v a d v a +u a d u a. (4.7) 

Rotations in the space of U(l) potentials (2.7) induces the corresponding 
transformations of electric and magnetic potentials: 

v a -> Vt a h v\ u a -> n a b u b . (4.8) 
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They are generated by 



K a b = v a d v b — v b d v a + u a d u b — u b d u a 



(4.9) 



S duality 

Four-dimensional S-duality (2.10) induces the corresponding transforma- 
tions of the target space variables. They consist of the axion shift 



k — > k + d 2 , u a — > u a + d 2 f a , 



proper electric-magnetic rotations 



u a ^u a , v a ^v a + u a di, z-^z^ + dx, 



and SL(2, R) scale transformation 



The corresponding Killing generators 

=d K + v a d u a , K d * = (e" 4 ^ - K 2 )d K + nd^ + u a d v a , 
K d3 =ds- 2nd K + v a d v a - u a d u a 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



satisfy sl(2, R) commutation relations 

~K d3 ,K dl ] =2K dl , \K d3 ,K d2 ] = -2K d \ \K d \K d A=K d \ (4.14) 

Ehlers— Harrison sector 

Non-trivial part of the isometry group includes 2p Harrison transformations 
whose action on real target space variables is rather involved (for p = 1 case 
see [0). Fortunately, more concise form exists in terms of complex variables, 
which we discuss in the next section. Here we consider only infinitesimal 
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transformations. Note, that they are not direct generalizations of those for 



p — 1 ||TT|. The first set (electric) is 



k h 1 = 2v a fd f + v a d 4 , + 2w a d K + (/e 2 * - v 2 ) d v a + 2v b (v a d v t + u a d ub ) 
+ {x~ u b v b + nfe 2 ^ d u a + (v a (x + v b w b ) + w a (fe 2<t> - v 2 )) d x , (4.15) 

where v 2 = (v a ) 2 , and we do not distinguish between upper and lower SO{p) 
indices. The main effect of these transformations is to generate electric po- 
tentials when acting on electrically uncharged configurations. Their magnetic 
counterpart reads 

k h 2 = 2u a fdf + {Kv a - w a )d^ + 2 (Kw a + v a e-^) d K + 
2u b (v a d v t + u a d u t) - (u b v b + X ~ Kfe 2(/} ) d v a + (fe~ 2(t> - u 2 + K 2 fe 2< ^) d u a 
+ (u a x + v a (u b w b - /e -2 *) + w a (nfe 2 ^ - v b u b )) d x , (4.16) 

where w a = u a — nv a , u 2 = (u a ) 2 . 

The last, Ehlers-type |22[ generator, which is produced through the com- 
mutation relations 

K?\K?*] =26 ab K E , (4.17) 

reads 

K E = 2f x d f + wWdtj, + (w 2 - v 2 e'^)d K 
+ (v b (x + v a w a ) + w b (fe 2 * - v 2 )) d vb 
+ (u b x + v b (w 2 + nv a w a - fe~ 2(t> ) + w b {nfe 2<t> - kv 2 - v a w a )^j d u b 

+ (x 2 ~f 2 + fv 2 e~ 2(t > + w 2 (fe 2(t > - v 2 ) + (v a w a ) 2 ) d x , (4.18) 

where w 2 = (w a ) 2 . Altogether these Killing vectors form a closed (p + 3)(p + 
4) /2-dimensional algebra. To identify it with SO(2,p + 2) let us consider a 
4 + p-dimensional space endowed with the following pseudoeuclidean metric: 

da 2 = G^dc = -(def - {de? + (der + ... + (^ +2 ) 2 , (4.19) 
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and denote the corresponding so(2,p + 2) generators as L^, where fi, v — 
0, 9, 1, ...,p+2 and a, b = 1, ...,p. Then the following correspondence between 
L^v and the above Killing vectors can be established: 

L ab = K ab , L 0a = ±(K^-K™), L 0p+1 = ±(K d *-K s ), 

L 0p+ 2 = \{K° - K d > - K d * - K E ), L ap+1 = l(K* + K?), 

L P+ i P+ 2 = \(K° - K d - + K d - + K E \ L ap+2 = 1{K* + K), 

L p+ie = ±(K0 + K d > + K d - - K E ), L ep+2 = - l -{K d * + K s ), 

L oe = ^{K 9 + K d - - K d > + K E ), L 0a = \{K^ - K e a ), (4.20) 

The target space (3.18) can now be identified with the coset manifold 
SO(2,p + 2)/(SO(2) x SO(2,p)). 

5 Holomorphic representation 

The coset manifold SO(2, p + 2)/(50(2) x SO(2,p)) may be parameterized 
by p + 2 complex coordinates and endowed with the Kahlerian metric jTT]] . 
These coordinates can be regarded as generalization of Ernst potentials of 
the Einstein-Maxwell theory. The set consists of the axidilaton z , p complex 
electromagnetic potentials 

<f> a = u a -zv a , (5.1) 

and the gravitational complex potential 

E = if - x + v a $ a - (5.2) 

In absence of vector fields (5.2) reduces to the original Ernst potential e = 
f + ix multiplied by i. However, when there are no scalar fields, E does not 
reduce to the Ernst potential of the Einstein-Maxwell theory. In particular, 
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i?-potential is not symmetric under an interchange v a <-> u a contrary to its 
Einstein-Maxwell counterpart. This reflects an intrinsic asymmetry of the 
N = 4 supergravity with respect to electric and magnetic sectors. 

In terms of complex variables the target space metric (3.18) can be rewrit- 
ten as 



dl 2 = Imz dE + 2Im$ a cM> a - (Im$ a ) 2 ^_ - i 



Imz 



Imz 



+ 



dz 



Imz 



where the quantity 



V = ImElmz + (Im$ a ) 
is related to the Kahler potential K as follows 



V 



-K 



The Kahlerian metric (5.3) is generated by K via 

dl 2 = 2 (d a dpk) dz a dz /3 , 
where holomorphic variables z a are enumerated as follows 
z a = (E } § a } z), a = 0,l,...,p + l. 



(5.3) 
(5.4) 

(5.5) 
(5.6) 
(5.7) 



Now the isometries of the target space may be presented as holomorphic 
maps z a — ► z a {z^) leaving K invariant up to an admissible transformation 



K -> K + h(z) +h{z), 



(5.8) 



where h is an arbitrary holomorphic function of the Kahler coordinates. First 
we describe two discrete maps which simplify substantially the derivation of 
holomorphic counterparts to continuous isometry transformations. One of 
them consists in the interchange of the Ernst potential E and the axidilaton: 



E 



z, z 



E, $ a unchanged. 



(5.9) 
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This transformation does not modify V (5.4) and hence leave the target space 
metric unchanged. In terms of the real variables it corresponds to 

f^jj^> e*-/-„V*, v^j^— 2 , W-x~« 

(5.10) 

where v 2 = {v a ) 2 , while u remains unchanged. The second discrete transfor- 
mation is more sophisticated. It corresponds to an interchange of the Ernst 
and axidilaton variables accompanied by the coordinate-dependent rescaling: 

E^h~ l z, z^hT x E, <$> a ^ h~ l ® a , (5.11) 

where 

h(z a ) = Ez + <$> 2 , $ 2 = ($ a ) 2 (5.12) 

is the function entering into the corresponding (admissible) transformation of 
the Kahler potential (5.8). This transformation is rather complicated being 
expressed in terms of real variables, but Kahler representation makes it easy 
to check. This map has another interpretation in the case p = 1 \T2, [T7] as 
inversion of the "matrix Ernst potential" . 

Now we are in a position to list holomorphic maps corresponding to con- 
tinuous isometries of the target space. Two simplest maps, linear shifts of E 
and $ a on real constants g and m a , correspond to gravitational 

E' = E + g, $' a = $ a , z ' = z, (5.13) 

and magnetic 

E' = E, $ /a = $ a + m a , z' = z (5.14) 

gauge transformation (4.3-4). An electric gauge transformation looks some- 
what more complicated: 

z' = z, $ /a = $ a + e a ^, E' = E — 2e a $ a — {e a ) 2 z, (5.15) 
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while the scale transformation (4.6) is simply rewritten as 

E' = e 2s E, $ /a = e s $ a , z' = z. (5.16) 

Now apply the discrete map (5.9) to the electric gauge transformation 
(4.2). Decomposing real Killing vectors of the previous section into holomor- 
phic and antiholomorphic parts one obtains 

K e a = zd^a -2<S> a d E , (5.17) 

so that after (5.9) we get 

K m =2® a d z - Ed^a, (5.18) 

what can be identified with the holomorphic part of the electric Harrison 
Killing vector (4.15). This means that the finite electric Harrison holomor- 
phic map may be constructed via an interchange of E and z in (5.15): 

E' = E, $ a = $ a + zAE z' = z - 2z/ a $ a - u 2 E, (5.19) 

where v a is a set of real parameters, and v 2 = (u a ) 2 . 

Similarly applying the discrete map (5.11) to the magnetic gauge (4.3), 
which infinitesimally reads 

K™ = d* a , (5.20) 

one obtains 

K^ 2 = (Ez + $ 2 )<9$ a - 2$ a (zd z + Ed E ) - 2$ a $% c . (5.21) 

This can be identified with the holomorphic part of the magnetic Harrison 
Killing vector (4.16), and thus applying (5.11) to the finite magnetic gauge 
transformation (5.14) one obtains the holomorphic map corresponding to 
finite magnetic Harrison transformations 

£' = f, $ /a = gl + ig A = l + 2fi a $ a + f i 2 (Ez + ® 2 ), (5.22) 
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where fi 2 = (/i a ) 2 . In the same way acting by (5.11) on the gravitational 
gauge 

K g = -d E , (5.23) 
one obtains Ehlers holomorphic generator 

K E = E (Ed E + $ a 9$a) - <$> 2 d z . (5.24) 

The corresponding finite form is 

E' = — — — , $' a = — , z > = z+ ^—. (5.25) 

1 + cE' 1 + cE' 1 + cE K ' 

The remaining ^-duality transformations can be obtained in a similar 
way by applying (5.9) to the gravitational gauge (5.13), scale (5.16) and 
Ehlers transformations. This results in 

z) E' = E, $' a = $ a , z ' = z + b, 

ii) E' = E + -^—, # = — _ , z>=-^—, (5.26) 
1 + d\z 1 + d\z 1 + d\z 

Hi) E' = E, = e d3 $ a , z' = e M3 2. 

Thus we have constructed [/-duality finite SO (2, 2 + p) transformations 
in terms of holomorphic maps acting on complex potentials which can be 
regarded as Ernst-type variables. In the case p — 1 our holomorphic maps 
coincide with given previously in where they were found in a different 
way using a symplectic representation of the so(2, 3) algebra. In the next 
section we will consider their application to solution generation purposes. 



6 50(2, 2+p) covariantization of the Schwarz- 
schild solution 

Target space of the vacuum Einstein gravity SL(2, R)/SO{2) can be param- 
etrized by a single complex variable E = if — \ and form a subspace of (5.3). 
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Using holomorphic transformations of the previous section one can construct 
a fully 5*0(2,2 + p) covariant counterpart to any vacuum solution of the 
Einstein equations. An application of this procedure to the Schwarzschild 
metric is likely to produce a general static black hole solution to dilaton- 
axion gravity with multiple vector fields. Brief results were reported recently 



23|| . Here we give more detailed derivation and discuss further properties of 



the solution obtained. 

In terms of holomorphic potentials the Schwarzschild solution reads 

E = if , /o = l-— , 2o = l, $S = 0, (6.1) 

and the three-metric (which remains non-affected by transformations) has 
non-zero components h rr = 1, h e9 = forf, h w = f r% sin 2 6. Let us first for- 
mulate the asymptotic conditions for holomorphic variables ensuring asymp- 
totically flat (Taub-NUT) behavior of the solution. It is worth noting that 
the NUT parameter enters as one of the SO (2, 2 +p) charges, so it has to be 
included in order to maintain the 5*0(2, 2 + p) covariance. 

If asymptotic values of the dilaton and axion are allowed to take arbitrary 
values 0(oo) = 0oo, k(oo) = k^, the charges should be defined as follows in 
the limit r — > oo: 

E^ih-m), (6.2) 
V r J r r 

where z^ = k m + ie~ 2 ^ x and three complex charges are introduced M. = 
M + iN, (ADM mass and NUT), V = D + iA, (dilaton and axion), Q a = 
Q a + iP a (electric and magnetic). Holomorphic transformations listed in the 
previous section can be combined into several sets preserving the asymptotic 
conditions (6.2) p3| . It is convenient to perform first the charging (Harrison) 
transformations and then Ehlers and one of 5-duality rotations imposing a 
condition = i, and afterwards restore an arbitrary value of via (5.26). 
So, as an initial step, we apply Harrison transformations (first magnetic and 
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then electric) combining them with suitable magnetic and electric gauge and 
axidilaton rescaling to preserve conditions E(oo) = z(oo) = i, $(oo) = 0. 
This results in 

imn f + 2q(l - /„) imn f - 2q(l - f ) 

k, = , z = , 

nm f nm f 

$a = (l-/o) (/* a n + »/»(2g-»m)) 
nrrif 

where 

q = v a fi a , n = 1 — z/ 2 , m = 1 — /i 2 , 
m f = 1- /j, 2 f , n f = l-u 2 f , rh f = f - fi 2 , h f = f -is 2 . (6.4) 

From here one can extract the following expressions for the metric function 
/ and the NUT potential: 

f = ^, x = ^^. (6-5) 
m f n f ° m f n f 

The latter has a non-vanishing Coulomb term in the asymptotic expansion, 
and hence, according to (6.2), the solutions is asymptotically Taub-NUT 
with N = 2qM /mn. (One can check that this is the only effect produced 
by X) short range terms in x are exactly compensated by the contribution of 
vector fields in (3.9), so the metric does not contain other terms than NUT 
due to u)idx l .) 

The drawback of this solution is that the NUT parameter is not free but 
is determined by charge parameters, and it is non-zero unless the SO(p) vec- 
tors fi a , v a are orthogonal, or one of them is zero. A NUT parameter of the 
solution, however, can be made arbitrary by applying the Ehlers transforma- 
tion (5.25), what can be easily seen considering its action on the asymptotics 
(6.2). Similar role in the S'-duality subgroup is played by the map (5.26 ii): 
it rotates dilaton and axion charges in the same way as Ehlers map rotates 
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the ADM mass and NUT. Combining both these rotations with gauge and 



scale transformations to preserve asymptotic conditions [23|, one obtains for 
the complex variable E the following expression 

E _ 2q(l - /p)(l + be) - n(cm } + brh f ) + im (h f - bcn f ) 
2g(l — /o)(c — b) + n(rrif — befhf) + im(brif + enf) 

where c is the parameter of the Ehlers transformation, and b = d.2 is that of 
(5.26 ii). The corresponding values of the ADM mass and NUT parameter 
depend only on c: 

M = ,f° 2 , {(1 - A 2 )(l - c 2 ) + Acq} , 
mn(l + c z ) L J 

iV = 2 ^ {c(l - /, 2 , 2 ) - (1 - c 2 )g} . (6.7) 

mn(l + c z ) L J 

Choosing c to be the root (regular at q = 0) of the equation 

g(c 2 - 1) + c(l - /i 2 z/ 2 ) = 0, (6.8) 



one selects strictly asymptotically flat solutions without NUT. In this case 

M= m±*) {l _„ v)i (6 . 9) 

mn(l — cr) v 7 

where 



i _ , 2,. 2 r 9 n 1/2 

c = -^^- + {l+((l-A*V)/2?) 2 } . (6.10) 

If M > 0, the ADM mass of the new solution is non-negative if /i 2 < 1, v 2 < 
1; c 2 < 1 (the limiting values should be considered more carefully in view of 
the singularity ). These restrictions on the parameters will be assumed so 
forth. Alternatively, by an appropriate choice of c one can construct massless 
solutions, these are non-trivial is N ^ 0. 

Finally, arbitrary asymptotic values of dilaton and axion can be generated 
via (5.25 i, iii), the metric being unaffected. The resulting complex axidilaton 
field reads 

z = K c -2^ M 1 - /o)( 1 + hc ) + njcfhf + bmf) - im(n f - bcn f ) , g ^ 
2q(l — /o)(c — b) + n(rrif — bcrhf) + im(brif + en/) 
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According to asymptotics (6.2) one obtains the following values of dilaton 
and axion charges in terms of the parameters introduced: 



mn(l + b 2 ) 
2M 



{b(S - /i 2 ) - (1 - b 2 )q} . (6.12) 



mn(l + b 2 ) 

Hence b plays for D and A the same role as c for ADM mass and NUT charge. 
In particular, any of scalar charges D, A may be set zero by an appropriate 
choice of b. 

The remaining set of U(l) complex potentials is given by 

$a = e-*~(l + c 2 y/ 2 (l + b 2 Y' 2 {l - f ) {fi a n + v\2q - im)} 
2q(l — fo)(c — b) + n(mf — bcrhf) + im(bnf + chf) ' 

while the corresponding electric and magnetic charges are 

V2M {(n/j, a + 2qu a )(b + c) + mv a (l - be)} 



Q a 



mn(l + 6 2 ) 1 /2(i + c 2)i/2 



_ V2M {(n^ + 2qu a )(l - be) - mu a (b + c)} 

mn(l + 6 2 ) 1 /2(i + c 2)i/2 • > 

Altogether the solution contain 2p + 5 arbitrary parameters M , /i a , u a , b, c, 
K oo, 0oo- However, the number of independent physical charges together with 
asymptotic values of dilaton and axion is 2p + 4. Indeed, both dilaton and 
axion charges (as can be directly checked using (6.7, 12, 14)) are determined 
by the equations 

M{(P a ) 2 - (Q a ) 2 } - 2NQ a P a 



D = 



2(M 2 + N 2 ) 



N{(Q a ) 2 -(P a ) 2 }-2MQ a P a 
A ~ 2(M 2 + N 2 ) ' (6 - 15) 

through electric and magnetic charges, mass and NUT parameter, or, in the 
complex form, 

(Q a ? 
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and thus are not independent parameters. Meanwhile the uniqueness theo- 
rem holds |24j] for harmonic maps to symmetric spaces (which is the present 
case too) saying that the solution is uniquely determined by the Coulomb 
charges. Moreover, the metric depends in fact only on three combinations of 
z/ a , fi a , namely /i 2 , u 2 , q = {i a v a (the only scalar invariants which can be built 
from SO (p)— vectors) . Also it can be shown that the metric is /^-independent. 
Extracting the function / from (6.6, 11, 13), one finds 

/ = h™<i + c 2 ) 

rUfUf + c 2 rifmf + 2cg(l — /q) 

that is the metric is determined by z/ 2 , /i 2 , q, c, Mq and a NUT parameter 
entering through \. All this indicates, that b is a gauge parameter. It is, 
however, useful to keep it arbitrary if one wishes to get solutions with a 
prescribed value of one of scalar charges. In particular, for 

' 1,1 +r^4|i + (V/(^ 2 -/i 2 )) 2 j (6.18) 



Z/ 2 — yU 2 |/i 2 — Z/ 2 | v J 

(as in (6.10) we choose the root which remains regular for all values of the 
parameters) the dilaton charge will be zero. 

The real potentials v a , u a , <f), k take simpler form in a "symmetric" gauge 
b = c. Then the scalar fields are 

e -2(^-<M = m i m f n f + c2 V^y + 2c g( x ~ ./o)l ( 6 19 ) 



n (m 2 + c 2 m 2 \ 



_ (/o ~ 1) {cm(/i 2 - u 2 )(l + f ) + 2q(m f - c 2 m f )} 

K — Koo — 7 ~ Z n\ ' (O.ZUJ 

n ym j + c z m j\ 



while the U(l) p potentials read 

a _ eM/o - !) {(c 2 m f - m f )u a - c(l + f ){^ a n + 2qu a )} 



2\ 



rrifUf + c 2 nfrrif + 2cq(l — f$ 



(6.21) 



w a = e^Hfo ~ 1) {(c 2 rh f - m f )(fj, a n + 2qv a ) + c(l + f )v a m 2 } 

n yrrif + c z mj\ 
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where we have chosen w a = u a — nv a instead of u a which are pure magnetic 
asymptotically when (p^ ^ 0, ^ 0. 

The solution (6.17, 19-22) constitutes the 5*0(2, 2 +p) covariantization of 
the Schwarzschild solution to which it reduces when all parameters /z 2 , u 2 , c, 
set zero. It contains a maximal number of [/-duality parameters 
compatible with asymptotic conditions and up to 2p trivial constants which 
can be added to v a , u a . It has arbitrary values of 2p vector charges, ADM 
mass, NUT parameter, and the scalar charges which are determined in terms 
of other charges via (6.15). This latter condition seems to be necessary for 
the absence of naked singularities (although no general proof has been given 
so far). Thus, by the uniqueness theorem, which holds for harmonic map into 
symmetric spaces [|4], it is likely to be the most general black hole solution 



(possibly endowed with NUT) to N = 4, D = 4 supergravity. 

7 Dilaton— Axion— Reissner— Nordstrom black 
holes 

Now let us discuss the choice of coordinates for the transformed metric. One 
reasonable choice is motivated by the Garfinkle-Horowitz-Strominger gauge 
for a (non-dyon) dilatonic black hole H, in which the curvature singularity 
sits at some finite radius r = r_: 

hoe = fori = fr(r - r_). (7.1) 

Taking into account that the three-space metric hij is not changed by trans- 
formations, one obtains for r an equation 

r(r -r_) = (r -rf)(r -Tq), (7.2) 

where 



M, 

ran 



^s 2a -0±^- iaj y /2 }, (7.3) 
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with 



a = c 2 -2cg + /iV, (3 = (/i 2 + u 2 )(l + c 2 ), 7 = 1 + 2cq + c 2 /rV. (7.4) 
Clearly 

r = r -ro, r_ = - rjT, (7.5) 
and hence the four-dimensional metric reads 

ds 2 = /(rft - 2iVdy?) 2 - /"Mr 2 - r(r - r_) (dfl 2 + sin 2 0<fy> 2 ) , (7.6) 

where 

/ = l-^ + ^ry> -K = rUro+2M). (7.7) 

Here M is a physical mass given by (6.7). The following relation between 
the seed and physical masses is useful: 

M -M=(r+ + r -)/2, (7.8) 

Unless K = 0, g 0Q is not Schwarzschildean, and the metric (6.28) does 
not coincide with that of ||. Hence in the general case the gauge (6.23) 
has no advantage, and one can choose more familiar curvature coordinates 
imposing a condition 

hee = for 2 = fp 2 . (7.9) 

Then the relation between the old and new coordinates is quadratic, and one 
obtains 

ro = I [4 + ro + (r 2 + 4p 2 ) 1/2 } . (7.10) 

In the curvature coordinates the transformed metric is rather similar to the 
Reissner-Nordstrom one 

& = (1 _ m. + °) {dt _ 2Niv? _ u _ m + or 1 dp , _ pW _ 

V p p J Vp p J 

(7.11) 
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with the same mass M and "charge" parameter 

C = r f(r - 2M ) + ^|(r + - 2M ). (7.12) 

It differs, however, from the usual Reissner-Nordstrom-NUT metric by a 
"variable mass" factor 

S = {l + (^/(2 P )) 2 } 1/2 , (7.13) 

which tends to unity when p — > oo. If r_ = 0, the solution is pure Reissner- 
Nordstrom-NUT. In this case, in view of (7.8), K = C, and hence (7.6) 
gives the same result. In what follows we will refer to the metric (7.11) as 
Dilaton-Axion-Reissner-Nordstrom-NUT (DARN-NUT) . 

Generally, the equation g 00 = has two solutions, p±. Taking into account 
(7.9), one can realize that smaller root corresponds to r = 0, while bigger 
one to r = 2M . This means that 

p 2 „=r -r+ p^ = (r --2M )(r + -2M ), (7.14) 

or, substituting (6.25), 

, _ 4M 2 (c 2 p 2 v 2 + l + 2cq\ 
P± ~ mn(l + c 2 ) I /iV + c 2 - 2cq ) ' {< Ab) 

Let us consider now the regular DARN black hole, assuming c to obey (6.10), 
so that iV = 0. Parameter r_ entering the above formulas then will be 

r_ _ {(1 - c 2 ) 2 (/i 2 - v 2 ) 2 + 4c 2 (l - p 2 v 2 ) 2 } 1/2 



2M (1 + c 2 )(l - n 2 u 2 ) 

while radii of the horizons are 



(7.16) 



2Mm 1 / 2 n 1 / 2 (l-c 2 ) 1 / 2 / (1 - c 2 p 2 u 2 f /2 \ 

P± ~ (1 + c 2 )(l - p 2 u 2 ) [ (p 2 u 2 - c 2 ) 1/2 ) ■ {7A7) 

The solution does not contain naked singularities if p 2 v 2 < 1, c 2 < p 2 v 2 
(what will be assumed below). 
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The inner horizon is pushed to the singularity when c 2 = fi 2 u 2 . In this 
case the solution is NUT-less if q = c. Both relations together imply that in 
this case the SO(p) vectors fi a , v a are parallel. It can also be checked that 
Tq = and hence K = 0, in which case we recover the standard expression 
for the metric of a "dilaton black hole" |§ with singularity at r = r_, 

2Mfc a + v 2 ) _ m 2 + (n 2 , 718 . 

l + /i 2 ^ 2 M ' 1 ] 

The electric and magnetic charge vectors are also parallel and in the gauge 
b = c read 

Qa V2M (l+^)u^ pa y/2M 0/ ^ ^ _ M (l + ^) 
nm m mn 

The dilaton and axion charges are still given by general formulas (6.15) (with 
N = 0) which now read 

(pay _ (Qa)2 M {( ^2 _ ^ _ ^2) _ ^2} 



D 



2M (l+/iV 
Q a P a _ 2M(1 + 

^ " " (1+^)2 ' ^ 2Uj 

where /i, ^ are norms of the vectors fi a , v a . Note, that in the case of a single 
(p = 1) vector field q = fiu inevitably, and the solution is NUT-less if q — c, 
what is just the present case. Therefore, the generalization of the "dilaton 
black hole" solution to the case of multiple vector fields is DARN solution 
with parallel SO(p) vectors /i a , z/ a , (or, equivalently, parallel Q a , P a ). It can 
be presented in the curvature coordinates as (6.33), where now 

C = Mr^ = (Q a f + (P a ) 2 , (7.21) 

while the radius of the event horizon is 

/ \ 1/2 
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Now let us clarify conditions, under which DARN solution becomes pure 
Reissner-Nordstrom. This happens when r_ = 0. From (7.16) it is clear 
that one should have \x = v, c = (thus implying q = for vanishing NUT). 
Therefore the only occurrence of the proper Reissner-Nordstrom solution 
in the dilaton-axion gravity with multiple vector fields is the case of the 
orthogonal SO(p) vectors /i a , v a of equal length. In this case 



1 — {l l V l 1 — {l l V l 



and the "charge" parameter in (6.11) assumes its standard form C = (Q a ) 2 + 
(P a ) 2 . From (6.19), (6.20) it is seen that in this case 

(f> = 0oo, K = Koo (7.24) 

("frozen moduli" [p5||). The fact that equal electric and magnetic charges be- 
longing to different U(l) sectors of N = 4 supergravity generate the Reissner- 
Nordstrom metric was discovered some time ago || within a truncated model 
with no axion. Our results show that this is the unique configuration of the 
N = 4, D = 4 supergravity (and its arbitrary-p generalization) ensuring such 
a property. In particular, if one takes c = b = (in which case we go back 
to the intermediate form (6.3) of the solution up to modifications due to an 
arbitrary asymptotic value Zoo), and q = to have N = 0, then 



1 — jJ, z V z 



and 

2M \ 2 /mV, /j 2 >u 2 
1-H 2 v 2 ) \n 2 /j, 2 , v 2 > ii 2 

so the metric is neither Reissner-Nordstrom, nor dilatonic if /i ^ v and 
/i, v 7^ 0. However, it is dilatonic if either all electric, or all magnetic charges 
are zero (then K = 0). From this formula one can see that there is another 
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case when K = 0, namely, when either fi or v is approaching unity. This is 
the BPS limit, which is worth to be discussed separately. 

Using formulas for physical charges in terms of the parameters one can 
show that the following identity holds generally 

M 2 + N 2 + D 2 + A 2 - (Q a ) 2 - (P a ) 2 = M 2 . (7.27) 

Thus the BPS limit of the solution corresponds to Mq — ► 0. From (6.9) it 
is clear that in the NUT-less case one possibility to achieve this limit (for 
non-zero M) is to take one of two quantities /z 2 , v 2 approaching unity and 
the other still keeping some different value. Then c ^ 1, and r_ — > 2M. In 
the limit one gets the following BPS saturation conditions 

(Q a ) 2 + (P a ) 2 = 2M 2 , D 2 + A 2 = M 2 . (7.28) 

(we choose the gauge b = c in (6.12)). From (7.17) it also follows that 
both horizon radii shrinks to the singularity. This is the "dilatonic" BPS 
state which corresponds in N = 4 supergravity (p = 6) to 1/2 unbroken 
supersymmetries 0. Indeed, absolute values of two central charges in our 
variables read 

M 2 = \ {(Q a ) 2 + (n 2 ± \/{Q a ) 2 (P a ) 2 - (Q a P a ) 2 ) ■ (7.29) 

For parallel charge vectors the square root vanishes, so both central charge 
moduli are equal, and from the saturation condition |^i )2 | = M with account 
for (7.20) one obtains (7.28). 

If both fi 2 , v 2 tends to unity simultaneously, then three cases should be 
distinguished: 

i) fi a , v a are parallel. 

In this case q = fiu, implying c = fiv to eliminate NUT, and we come 
back to the general dilatonic case discussed above. 

ii) /j, a , v a are neither parallel, nor orthogonal. 
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From (6.10) one finds that then c — > 1 so that 

^^-g, g^o, (7.30) 

and then r_ — > 2M, p± — * 0, hence the limit is also dilatonic. 
hi) fi a , v a are orthogonal. 

In this case q = 0, c = 0, and both horizons merge at 

2My/mn 

P± = Pext = ~, 2~2- (7.31) 

The limit of the right hand side depends on the curve in the /x, v plane along 
which the point /x = 1, v — 1 is reached. If one takes the limit along the line 
\i = v, one gets 

Pext -> M, (7.32) 

what corresponds to the extreme Reissner-Nordstrom black hole. If one goes 
to the limiting point at some different angle, e.g. 

(j, — 1 — e, v — 1 — ke, e — > 0, (7.33) 

one generally gets an extremal solution with the horizon radius 

PL = TTT^- (7-34) 

In two limiting cases, k = 0, k = oo (non-dyons) the radius tends to 
zero, hence solution is dilatonic. For k — 1 it is the extremal Reissner- 
Nordstrom, for all other values of k this is an extremal limit of the generic 
solution (7.11). Indeed, the parameter r_ determining the deviation from 
the Reissner-Nordstrom geometry now is given by 

the charge parameter is C = (Q a ) 2 + (P a ) 2 , while its ratio to the square of 
mass is 

{Qa) 2 ± { pay 2(1 + fc2) 
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The dilaton charge is equal to 

(P") 2 -(Q") 2 _ 2M 2 (fc 2 -l) 

while the axion charge is zero, so the force balance condition is satisfied 

M 2 + D 2 = (P a ) 2 + (Q a ) 2 . (7.38) 
Furthermore, the electric and magnetic charge vectors have the form 

where now /i a , v a are unit vectors. Hence k can be expressed through the 
norms of charge vectors Q = (Q a Q a ) 1/2 , P = (P Q P a ) 1 /2 as 

* = ^ (7.40) 

From (7.38) it also follows that 

P + Q = V2M, (7.41) 

and then, in view of (7.36-37), 

P-Q = V2D, (7.42) 

thus only one BPS bound of iV = 4 supergravity is saturated, signalizing 
that the solution has N — 1 residual supersymmetry ||, |J. 

Therefore, in the BPS limit the DARN solution with orthogonal charge 
vectors has the horizon radius (7.33) taking all values in the interval 

< p ext < M. (7.43) 

Zero value is reached for purely electric k = or purely magnetic k = oo 
configurations, while the upper limit corresponds to a symmetric dyon Q = 
P. 
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To avoid confusion, note that this interpretation is based on our choice of 
curvature coordinates (like in the standard Reissner-Nordstrom metric). In 
[P] it was claimed that all BPS dyonic black holes in purely dilaton gravity 
with two vector fields, possessing electric and magnetic charges in different 
U(l) sectors, have the same radius of the horizon equal to that of a pure 
Reissner-Nordstrom black hole. This may seem to disagree with our conclu- 
sion since the particular model of f| lies within the scope of our consideration. 
There is no contradiction, however, since in [f| a different coordinate system 
was used, namely their radial variable (denoted by prime) is related to our r 
(7.1) as r' = r — r_/2. Taking into account (7.35), one can easily show that 
in the BPS limit r' = M independently on k. But the geometrical radius of 
2-spheres in terms of primed variable is equal to {r' 2 — (r_/2) 2 } 1 / 2 , whereas 
it is p 2 in our case. Thus one has the same expression for the area of the 
horizon surface: 

A H = 47T/4, = SvrPQ. (7.44) 

Recall that this expression is valid only in the case of orthogonal charge 
vectors. If desired, the solution (7.11) may be presented in the form similar 
to |J. However, we prefer to use curvature coordinates because of their more 
transparent geometrical meaning. 

8 Conclusions 

We have given a concise representation of the three-dimensional [/-duality 
for the (generalized) bosonic sector of the stationary N = A supergrav- 
ity in terms of Kahler parameterization of the target manifold SO(2,p + 
2) /(SO (2) x SO(2,p)), where p is a number (chosen arbitrary for gener- 
ality) of U(l) vector fields. Global 5*0(2, p + 2) symmetry of the three- 
dimensional theory was described in terms of holomorphic maps exhibiting 
the target space isometries. They include Ehlers-Harrison transformations, 
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S-duality, as well as SO(p) rotations, gauge and scale transformations. Se- 
lecting transformations which preserve an asymptotic flatness, one obtains 
a (2p + 4)-parametric subset which is suitable to perform a SO(2,p + 2) 
covariantization of any asymptotically flat solution to the vacuum Einstein 
equations. Such covariantization of the Schwarzschild solution presents the 
most general static black hole in the dilaton-axion gravity with multiple 
vector fields. 

Generically DARN black hole solution obtained possesses two horizons. 
The metric depends on four parameters (in the NUT-less case): the ADM 
mass, and three SO(p) invariants built out of two SO(p) charge vectors. For 
parallel charges an internal horizon shrinks to the singularity and the metric 
coincides with that of the standard dilaton black hole. For orthogonal charge 
vectors of equal length one has a pure Reissner-Nordstrom metric. All other 
charge configurations correspond to generic DARN metric with two non- 
singular horizons and a space-like singularity. The NUT-generalization of 
the solution is also given. 

Having now general non-extremal solution, one can better understand 
how generic are particular configurations found earlier in the Bogomol'nyi- 
Gibbons-Hull limit |S], [7|, |], One finds that in the BPS limit general 
DARN solution has a "dilatonic" behavior (both horizons shrinking to the 
singularity), unless charge vectors are orthogonal. Thus the generic BPS 
static black hole in the N = 4, D = 4 supergravity has a vanishing area 
of the horizon (and N=2 residual SUSY). In the case of orthogonal charge 
vectors two horizons merge at some (generically finite) radius depending on 
the ratio of norms of the electric and magnetic charge vectors. Non-dyon 
configurations are still dilatonic, while dyons have finite radius of the horizon, 
taking maximal value in the "symmetric" case. 

Few remarks are in order concerning the relationship between our solution 
generating technique and those, frequently used in the string theory, which 
are usually based on the presentation of the four-dimensional theory as di- 
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mensional reduction of some higher-dimensional one. First of all, our direct 
construction of the sigma-model does not require the theory to be obtainable 
via dimensional reduction. Also, the standard reduction procedure usually 
gives rise to some matrix representation in terms of real variables which is 
manifestly invariant under action of the symmetry group, but physical in- 
terpretation of transformations often remains obscure. The advantage of the 
present approach is that it deals with complex potentials, thus reducing the 
total number of variables by half. It also provides from the very beginning 
for a clear physical identification of transformations involved. In addition, 
as we have shown, some discrete symmetries exist which are manifest only 
in terms of Kahlerian variables. 
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